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The exact form factor for the 0(3) non-linear a model is used to predict the three-magnon 
contribution to the spin correlation function, S(q,u>), near wavevector q = tv in an integer spin, 
one-dimensional antiferromagnet. The 3-magnon contribution is extremely broad and extremely 
weak; the integrated intensity is < 2% of the single-magnon contribution. 

PACS number 75.10.Jm 
The Hamiltonian for the one dimensional Heisenberg antiferromagnet of spin s is 

H = jJ2Si-S i+1 . (1) 



i Based upon the large-s limit we write the spin operators as 



Sj««(-i) j >Cj) + *Cj). (2) 

where s(f> and I represent the staggered and uniform magnetization of the spin chain. We set the lattice spacing to 1 . 
The low energy behaviour of this Hamiltonian can be described by the 0(3) non-linear a model [[[). Recently, the exact 
form factors for this field theory were calculated [|]-§}. The resulting prediction for S(q,cu) at q sa was discussed 
in ||. Here we comment on the prediction at q « tt. In the continuum approximation, the (zero temperature) spin 
Q\ ' correlation function at<7 = 7r + fc«7ris 

on : 
■ 



S ab {n + k,uj) = s 2 dx / dtexpi(ujt - kx) < CI \ 4> a (x,t)4> b (0,0) \ >= 8 ab S(-K + k,ui), (3) 



where |fi > is the groundstate. 



We insert a complete set of asymptotic states. It is known that the spectrum consists of a triplet of massive 

o 
o 



magnons. Thus, the asymptotic states will be characterized by a spin index and a momentum for each particle in the 
state. 

| n >=| ai,pi;a 2 ,p 2 ; ■ ■ ■ ;«n,Pn > 

It is convenient to label the particles' momenta by the rapidities, 9i~. 

Ei = A cosh 0i, pi = (A/v) sinh 0, (4) 

where v is the spin- wave velocity corresponding to the velocity of light in the quantum field theory and A is the gap 
corresponding to the rest mass energy in the quantum field theory, (v w 2.49J and A w .4107J for the s=l chain.) 
Thus we may write: 

S ab {-K + k,w) = s 2 I dx I dtexpi{ujt- kx)J2 — T[ / T~ < n I ^ a ^ x ^) I nxn | fc (O,O) | ft > . (5) 

J J n U - i=l J 4?F 

We use 

< n | <j){x,t) | n >=< Q | 0(0) | n > exp[-i(E n t - P n x)} , (6) 
where P n and E n refer to the total momentum and energy of the state \n >, to obtain 

S(7r + fc,c) = ,s 2 ^££lf[ J ^5(k -P n )5(u,-E n )\<n\^ (0,0) \ n >\ 2 . (7) 
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The field is renormalized 



(8) 



in order that we satisfy the relation < fl | $°(0) | b,p >= S ab . Symmetry arguments guarantee that only asymptotic 
states with an odd number of magnons will offer non-zero matrix elements. 
For the one particle contribution, 



Si (it + k, ui) = s vZtt 



, „ S(uj - \/v 2 k 2 + A 2 ) 



y/v 2 k 2 + A 2 



The integrated intensity is, 



Si(q) = / dwSi(q,w) — s 2 vZir 



1 



\/v 2 k 2 + A 2 



(9) 



(10) 



Numerical simulations on the s = 1 antiferromagnet || indicate that Z « 1.26. 
The three particle contribution can be written: 



£3(77 + k,uj) = s 2 vZ 

1 

3 



d(9i - 9 2 )d(0 2 - 63) /—z — now 
-—2 6(\/uj 2 - v 2 k 2 ~ M(9i, 2 , 03)) 



y/uj 2 — v 2 k 2 Jo 
i |< I * a (0,0) I ai,pi;a 2 ,p 2 ;a 3 ,p 3 >| 2 , 



(11) 



a,ai ,a2 ,03 



where: 



M 



^ Pi and 



(12) 



»=i 



Remarkably, the 3-particle form factor has been calculated exactly using the integrability of the non-linear a model 



SI: 



I ]T |<r!|$ a (0,0)|a 1 ,p 1 ;a 2 ,p 2 ;a 3 ,p 3 >| 2 = tt 6 ^, 2 , 3 )\ 2 [2(0 2 2 i + 0\ 2 + ^1) + ^A, 

i /j(0i,9 2 ,e 3 ) = l[i>(0 ij ), 



a,ai ,a2 ,0,3 



i>j 



0(2ni 



ITT 2 

tanh — . 

-9) 2 



(13) 



The resulting integral in Eq. (11) can be easily performed numerically. The result for S(n,uj) is shown in Fig. 1. 
A S(7i,ra) 
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FIG. 1. 1-magnon and 3-magnon contribution to spin-correlation function, S(n,u>). The factor of s 2 vZit/A in Eq. (§) has 
been divided out so that the peak at uj = A has unit integral. 
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The three-magnon contribution vanishes below 3A. In the limit to — > 3A, it behaves as: 

63(77, u) -> s 2 vZ x .01045(w - 3A) 3 /A 5 . (14) 

It has a rounded, asymmetric peak atw« 6.33A then decays at high energy as s 2 vZ x 19.9/{w 2 [ln(w/A)] 2 }. The 
integrated intensity of the three particle contribution at q — ir is 

63(77) w .01935i(tt). (15) 

This 3-particle contribution to S(q, lo) is very weak and very broad. It is instructive to calculate the average frequency 
of the 3-particle term: 

- _ J dujujS 3 {n,uj) 
J do>6 3 (77,o;) 

Using the result of Q for the integral in the numerator of Eq. ([l6]) we find: 

o> 3 w 75. 2A. (17) 

This enormous value is explained, to some extent, by the rapid decrease of 63(77,0;) at to — ► 3 A and its slow drop 
off at large lo, making the integral in the numerator of Eq. (|l^) only logarithmically convergent. The contributions 
to 6 from still higher numbers of particles have also been considered . Indications are that these are still more 
negligible than this tiny 3-particle term. 

It is clear that this result for 53 cannot be applied completely to the s=l antiferromagnetic chain. In particular, 
from the numerically determined single magnon dispersion relation, which is quite well fit by: 



E(p) w -y/A 2 +w 2 sin 2 p, (18) 

it can be seen that the maximum possible energy of a 3 magnon state with total crystal momentum 77 is only about 
17A. 6(77,0)3) in the field theory gets significant contributions from bosons with momenta considerably larger than 
77, the maximum possible in the lattice model. In fact the relativistic approximation to the dispersion relation seems 
to break down significantly for p > .277. We note that, for u) < 9A, only magnons with p < .2n contribute to S3 (it, lo). 
Thus we might hope that S 3 (q, ui) calculated from the field theory is fairly accurate for 77 — q < .2tt and uj < 9 A. The 
non-relativisitic corrections to the magnon dispersion relation make 6*3(77,0;) vanish for lo greater than about 17A. If 
we only integrate over 83(17, ui) up to lo = 17A this reduces 53(77) to about .0125*1(77). However, the non-relativistic 
corrections may also tend to increase 63(77,0;) for 9 A < lo < 17 A since they flatten the dispersion relation hence 
increasing the density of states. Clearly the value of C03 in the s=l chain must be less than the maximum possible 
frequency of about 17A. If we make a rough estimate that 0)3 w 10A and 63(77) w .02S*i(77) then we can estimate the 
overall average frequency (also ignoring 5 and more magnon contributions) as: 

f dw\Si(iv, lo) + 63(77, lo)]lo . 63(77) . . . . . 

^ J r , L r'V TT^T 1 «A + -^a> 3 «Al + .02xlO =1.2A. (19) 
J oto;[6i(77, lo) + 5 3 (77, lo)\ 61(77) 

This estimate can be compared to the numerically determined value. We may use an exact sum rule for the Heisenberg 
antiferromagnet : 

' dLOLoS{q,Lo) = -(1/2) < [[H,S z (q)],S*(-q)] >= (2/3) |e | (1 - cos q), (20) 

where eo is the groundstate energy per site. Using the numerically determined value of S(q) and eo gives Lo(q). This 
is plotted in Fig. (18) of Ref. || where it is referred to as wgjyfj^. From this figure we see that at q — 77, Co ps 1.2A 
in agreement with our crude estimate of Eq. (|l^). This lends some confidence to our prediction that the 3-magnon 
contribution to 5(g, lo) near q w 77 is extremely broad and extremely weak, with 0)3 of order the maximum possible 
3-magnon energy, w 10A and the relative integrated intensity of order 2%. As such, it will be extremely difficult 
to observe experimentally. We note that the 2-particle contribution to S(q,Lo) near q — is also extremely difficult 
to observe |l]||f7| since it appears to only occur for q < .277 and since S(q,Lo) oc q 2 as q — > 0. Thus experimental 
applications of the beautiful exact results ||0 Q on the non-linear a model remain elusive. 

If we assume that the 3-particle form factor in Eq. (^) goes to a non-zero constant at the top of the 3-particle 
continuum, lo w 17A, then it follows from phase space considerations that 63(77,0;) drops discontinously to at that 
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energy. This is quite different than the continous vanishing of S3 at the bottom of the 3-particle continum which is 
entirely due to the vanishing of the form factor there. Possibly this discontinous drop at the top of the continuum 
might be easier to observe experimentally than other features. However clearly the tiny size of the drop would make 
even this extremely difficult. 

We would like to thank Bill Buyers for interesting us in this problem. This research was supported in part by 
NSERC of Canada. 
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